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We propose a scheme with dc-control of finite bandwidth to implement two-qubit gate for super- 
conducting flux qubits at the optimal point. We provide a detailed non-perturbative analysis on the 
dynamic evolution of the qubits interacting with a common quantum bus. An effective qubit-qubit 
coupling is induced while decoupling the quantum bus with proposed pulse sequences. The two-qubit 
gate is insensitive to the initial state of the quantum bus and applicable to non-perturbative cou- 
pling regime which enables rapid two-qubit operation. This scheme can be scaled up to multi-qubit 
coupling. 
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INTRODUCTION 



Superconducting Josephson junction (J J) qubits (for a 
review, see e.g. 

[I i, S II) provide an arena to study 
macroscopic quantum phenomena and act as promis- 
ing candidates towards quantum information process- 
ing. For the three basic types of superconducting qubit, 
namely charge qubit, flux qubit and phase qubit, sin- 
gle qubit coherent operations with high quality factor 
have been demonstrated in many laboratories @, 0, 
H, B 113, im, Il3l- However, the best way to achieve 
controllable coupling and universal two-qubit gate are 
still open questions. A number of experimental at- 
tempts [H M, m, M, H Eli, 19 , 2Q, 11, 2Z, 23., M 0, 0, 




as well as theoretical proposals have been put for- 

H lU, M, m, S, M, m, M, M, M, H, EH, Sii 

according to the characteristics of each specific circuit. 
In this paper, our discussion will be focused on coupling 
superconducting flux qubits 43, 3, 3]. 



The straightforward consideration to realize two-qubit 
entanglement is utilizing the flxed inductive coupling be- 
tween two flux qubits. With tunable single-qubit energy 
spacing, this fixed coupling can be used to demonstrate 
two-qubit logic gate [26| . However, tunable coupling is 
required to achieve universal quantum computing. At 
early stage, dc-pulse control is widely adopted in the 
tunable coupling proposals Hull. Main disadvantage 
for this method is the inefficiency to work at the degen- 
eracy point which is a low-decoherence sweet spot. At 
the optimal point, the natural inductive coupling is off- 
diagonal in the diagonal representation of the free Hamil- 
tonian. Hence the coupling only has second-order effect 
on the qubit dynamic for the detuned qubits. Another 
difficulty related with dc control is the operation error re- 
lated to the finite rising-and-falling time of the dc-pulse. 
Recently, more attention is paid to coupling schemes with 
ac-pulse control [H, H M, S 113, So] . While most of 
the ac-control coupling schemes can work at the degener- 



acy point and no additional circuitr y is needed 3J, |40| , 
some of them require strong driving |34l | or result in slow 
operation [i^. Meanwhile, unwanted crosstalk is present 
due to a always-on coupling. The possible solution to 
the above problems is the parametric coupling scheme 
with a tunable circuit acting as coupler (36| . A third flux 
qubit has been demonstrated as a candidate for this cou- 
pler [25, 37 1 . However incorporating additional nonlinear 
component to the circuit would increase the complexity 
of the circuit and might introduce additional noise and 
operation errors. 

In this paper, we propose a scalable coupling mech- 
anism of flux qubits with four Josephson junctions in 
two loops (4JJ-2L). The coupling is induced by a com- 
mon quantum bus, such LC resonator or a one- 
dimensional superconducting transmission line resonator 
(TLR). The effective coupling Hamiltonian is diagonal 
with the free Hamiltonian of single qubit at the optimal 
point. With appropriate dc-control pulse, a dynamic two- 
qubit quantum gate can be realized for superconducting 
flux qubits at the optimal point. The on-and-off of the 
coupling can be switched by dc-pulse of flnite bandwidth 
without introducing additional error. This protocol is 
based on the time evolution of a non-perturbative inter- 
action Hamiltonian. Therefore it is applicable to "ul- 
tra strong coupling" regime, where the coupling strength 
is comparable to qubit free Hamiltonian. Contrarily to 
parametric coupling which requires a strongly non-linear 
coupling element, the scheme described in this paper uti- 
lizes a linear element. The strong non-linearity of para- 
metric coupler required in order to achieve fast enough 
two-qubit gates induces strong imperfections of the gates 
and added the difficulties related with microwave con- 
trol [s^l- While the two-qubit gate based on linear cou- 
pler is intrinsically free of errors if proper DC control is 
achieved. Thus the linearity of the coupler considered in 
this work has not only the advantage to be insensitive 
to the state of the coupler, but also offers the possibil- 
ity of error- free gates. Due to these advantages, this new 
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proposal could be a promising alternative in experiments. 

This paper is organized as follows. In Sec. II, we first 
analyze the e nerg y spectrum of the 4JJ-2L qubit con- 
figuration 43|, |4J]. In Sec. Ill, the setup of our coupling 
mechanism for this type of qubit is described and the sys- 
tem Hamiltonian is derived. In Sec. IV, we present two 
different pulse sequences to realize the effective two-qubit 
coupling and construct two-qubit logic gates. The char- 
acteristics of this coupling scheme based on experimental 
consideration are analyzed in Sec. V. The discussions of 
this paper are given in Sec. VI. 



circuit: 



(i) 
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where $m the magnetic flux threading the main qubit 
loop. Adding up the two equations in ([2]), we get 



(0 



(3) 



FLUX QUBIT WITH TUNABLE QUBIT GAP 

A single flux qubit discussed in this paper is shown in 
Fig. (HI). Each qubit is composed of four Josephson junc- 
tions in two loops: the main loop (lower loop) and the dc 
SQUID loop (upper loop). The main loop encloses three 
junctions: two identical junctions with Josephson energy 
Ef and one shared with the dc SQUID loop with Joseph- 



son energy a^'^ E^f' where a)^' is the ratio of the Joseph- 
son energy between the first two junctions and the third 
one (here and hereafter, the superscript (i) denotes the 
variables of the z-th qubit). The main loop forms a fiux 
qubit whose energy eigenstates are the superpositions of 
the clockwise and the counterclockwise persistent current 
states [43, 44 1- The 4- J J flux qubit is different from the 
conventional design of a fiux qubit due to the additional 
dc SQUID loop. The third junction of 3-JJ flux qubit 
is replaced by a dc SQUID in this 4-JJ design. There- 
fore the effective Josephson energy of the third junction 
can be controlled by the magnetic flux threading the 
dc SQUID loop. Assuming the two junctions in the dc 
SQUID loop are identical, the effective Josephson energy 



is a(^)($^'^)^f = 



2q!o cos 



E^^ with $0 the 



flux quantum. This feature, as we show later, enables the 
qubit gap to be tunable. This increases the in situ con- 
trollability of the quantum circuit 4^, 44 1 . The main loop 
and the dc SQUID loop of each qubit can be controlled 
by external on-site flux bias separately. A high-fldelity 
two-qubit operation has been proposed recently for the 
4JJ-2L qubit 

As shown in Fig. H]), the Josephson phase differ- 
ences of the four junctions in one qubit are denoted 
by ^p^i\ 'fi^ ) 'fi'^ and (^4^ respectively. By deflning 



f3 = 



in ' + "PI 

one qubit loop is 



)/2, the total Josephson energy in 



($|j''')i?j'-' cos (^3''', where we have used the 
quantization relation in the dc SQUID loop 



E 



(-0 
cos kp\ 



E 



COS + 

fluxoid 



(1) 



There are two other fluxoid quantization relations for this 



where ^[^^ = + ^'"^^ /2 is the total magnetic flux 
threading the qubit loop. Then the total Josephson en- 
ergy of the four junctions in the loop is 



a«($«)i?f COS (^2.^-(^«+^^^^)j 
+Efcosip'^;^+E^^cos^i'\ (4) 



It takes the same form as that of the 3-JJ flux qubit [43j, 
441 except that the ratio a ') is tunable. If the total 
magnetic flux is close to half a flux quantum $0/2 
and a^*^ > 0.5, the function Uq (v'i\^2^^ represents a 
landscape with periodic double- well potentials. 

With external flux bias, one can set the operation point 
in one double-well potential. The classical stable states 
of this potential correspond to the clockwise and the 
counter-clockwise persistent current states. By chang- 
ing the ratio a*-*' between the Josephson energy of the 
third junction (through the dc SQUID), the height of the 




FIG. 1: (Color online) The schematic of a single qubit with 
four Josephson junctions (denoted by cross) connected in two 
superconducting loops. The upper loop forms a dc SQUID 
with two identical junctions while the lower loop encloses 
three junctions similar as the conventional 3-JJ flux qubit. 
Each loop can be controlled separately by external magnetic 



flux $j ' and "J> 



respectively. 
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tunneling barrier (hence the tunneling rate) between the 
two minima of each double-well is tunable. When a^'^ is 
set in appropriate range, coherent tunneling between the 
two wells of the potential is enabled while the tunneling 
between different potentials is highly suppressed. 

Taking into account the electric energy stored in the 
four capacitors, we can get the full Hamiltonian of this 
system. The energy spectrum of the circuit with a*-*-* = 
0.8 and Ef /E^^'> = 35 {e'^ = /2C denotes the 
Coulomb energy of the first (second) junction of the i- 
th qubit and C is the junction capacitance) is shown in 
Fig. [5] as a function of the rescaled total magnetic flux 
/(») = $['V*o- In the vicinity of «>^*^ = $o/2, the lowest 
two energy levels are far away from other energy lev- 
els and form a two-level subspace which can be used as 
a flux qubit. The eigenstates of the flux qubit are su- 
perpositions of the clockwise and the counter-clockwise 
persistent current states. The 4-JJ flux qubit works the 
same as its 3-JJ prototype except that the barrier height 
of the double-well potential is tunable in situ. In the two- 
level subspace, the free Hamiltonian for the i-th qubit is 
written as 



7. <^;. ' H t:. ' 
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FIG. 3: (Color online) The energy gap of single qubit A^'^ 
and its derivative dA'^'/da''-' as a function of a'*' for (a) 



(5) 



Ej/E^^' = 35 and (b) Ey'/E)^' = 50 (for simplicity, the 
superscript (i) is omitted in the figure). The solid line (red) 
is obtained from the exact diagonalization of the original 4- 
JJ qubit Hamiltonian while the dashed line (black) is ob- 
tained from the analytical solution of the tight-binding model 
with WKB approximation which breaks down at low barrier 
regime. The energy is in the unit of Ej. 



where e^^'^ is the energy spacing of the two classical cur- 
rent states 



p 



$0 

2 



(6) 



and A*^') is the energy gap between the two states at the 



degeneracy point ^[^'' ~ ^o/2. 



AW($«) 



AW(a« 



2a. 



cos(^-^)). 

<J>0 



(7) 



can 
3] as 

W>*(i) _ 



According to the tight-binding model, A^'^ 
be evaluated through WKB approximation 
A(^) « (wa/27r)exp(-[4a(l + 2a)£;f /i;^'']i/2(s 
(/)*^'V2a)) where Ua is the attempt frequency of escape 
in the potential well and cos^*'*-' = O.Sa*-'^ (the Planck 
constant h is set to be 1). In Fig. [31 the energy gap A^*' 
and its derivative dA*^*' /da'*' arc shown as a function of 
a^^\ The results are obtained from numerical calculation 
and analytical derivation based on WKB approximation. 
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FIG. 2: (Color online) The energy spectrum of the lowest six 
energy levels of the superconducting loop with respect to the 
total magnetic flux . The energy is in the unit of Ej while 
the magnetic flux is in the unit of $o. We take = 35 

and Q<'' = 0.8. 



THE COUPLED SYSTEM 

A schematic to illustrate our coupling mechanism is 
shown in Fig. [4] with two different types of data bus, i.e., 
LC resonator and ID TLR. For simplicity, we first con- 
centrate on coupling two qubits. The problem of scale- up 
will be discussed later. As we described in the previous 
section, each qubit is composed of four Josephson jimc- 
tions in two loops: the main loop (the lower loop) and the 
dc SQUID loop (the upper one). The main loop and the 
dc SQUID loop of each qubit can be controlled by exter- 
nal on-site flux bias independently. The two qubits are 
placed in sufficient distance so that the direct coupling 
can be effectively neglected [IHj . The two qubits are both 
coupled with a common data bus such as a twisted LC 
resonator or ID on-the-top TLR via mutual inductance. 

Due to the mutual inductance with the resonator. 
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(a) 



QB2 



by the quantum bus in the upper loop and the lower loop 
are opposite. Inserting Eq. pO|) into Eq. ([8]), we find the 
total flux is contributed only by the external applied 



magnetic flux as $J: ~ <!>, 



,(0 . 

q,o 



Since the a, com- 



ponent of the qubit is coupled with the resonator 
contributes a pure ax coupling with no CTj component. 
Therefore the qubit can always be biased at the optimal 
point <^['^ = $o/2. 

For the quantized mode of the resonator, 



/=/^(a + at). 



(12) 



(b) 



FIG. 4: (Color online) The circuit design examples to im- 
plement the required coupling. Two 4JJ-2L flux qubits are 
coupled with each other through the inductive coupling with 
a resonator as data bus: (a) a twisted LC resonator and (b) a 
ID superconducting transmission line resonator in a separate 
layer. The current of the data bus induces magnetic fluxes 
both in the upper loop and in the lower loop of each qubit. 
The directions of two magnetic fluxes are opposite. 



where uj = {LC)~^^^ is the plasma frequency of res- 
onator, L (C) the lumped or distributed inductance (ca- 
pacitance) of the resonator and (a) the plasmon cre- 
ation (annihilation) operator. With these denotations. 



(i) 



where 



Jd 



M. 



-at) 



(13) 



(14) 



the magnetic fluxes include the contribution both from 

(i) 

the external applied flux and the resonator, i.e., 'I'm = 



1>Sb and 



o + b' where the subscript 
e(6) indicates the contribution from the external mag- 
netic flux (the quantum data bus) respectively. Then 
the total magnetic flux ^[^"^ reads 



(0 _ 



(0 



m,b 



d.b 



(8) 



The coupling between a single qubit and the data bus 
includes two parts: the coupling of the qubit with the 
dc SQUID loop via mutual inductance M^'"* and the cou- 
pling of the qubit with the main loop via Mm"*. The 
magnetic flux induced in the dc SQUID loop and the 
qubit main loop are 



d,b 



(9) 



respectively and / is the current in the resonator. For 
our purpose, the two magnetic fluxes satisfy 



d,b 



-2$ 



(0 



,b- 



(10) 



This can be implemented by designing the mutual induc- 
tance 



M 



-2/Vf(*) 



(11) 



The minus in ([TO)) is due to the special layout of the data 
bus so that the directions of the magnetic flux induced 



Usually the mutual inductance of the resonator and the 
qubit loop is about several pH to several tens of pH. 
For example, if we take M^*-* = 10 pH, w = 1 GHz and 

L = 100 pH, /^'V^o « 5.6 X 10"'' < 1. This means the 
magnetic flux contributed from the resonator is much 
smaller than that from the external applied magnetic 
field. To the first order, the energy gap of a single qubit 
is modified by the resonator as 



A ( a« 



A 



dA (a(*)l 



(5a«(a-Hat) 
(15) 



with 



aW = 2a« cos In^j 
fcW = 2a^*)^sin U^ 



(16) 



The Hamiltonian for a single qubit linearly interacting 
with the data bus reads, 



($«e + <i /2) A«(<f>g 
1^ . ^rr(') J 1^ L 



with 



(17) 



(18) 
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and the coupling coefficient 



with 



(19) 



(20) 



Note that magnitude of the coupling g^'^^ increases with 



the mutual inductance 



^dC). If $^^,e + *^ye/2 = (n + 

0.5)$o (where n = 0, ±1,±2 is an arbitrary integer), 
qubit is biased at the degeneracy point and the system 
Hamiltonian is written as 



(i) 



H 



ioa'a 




A«($«) 



By tuning the external magnetic flux ^ in the dc 

SQUID loop to be n<I>o, 5^*^ i^d c) = 0' t^^^ qubit is decou- 
pled from the resonator in the first order. The qubits act 
independently and single-qubit operation can be imple- 
mented by biasing <I>q*e together with microwave pulse. 

In the above discussion, the condition Eq. (fTTj) is as- 
sumed. However it might not be precisely satisfied in 
practical case. Suppose there is a small deviation in the 
fabrication process that = -2(1 + S)Mil^ (where 

6 <^ 1), the total magnetic flux ^[^^ includes a small 
contribution from the resonator, 




t q,c 



(0 

d,o 



(22) 



series of specific unitary evolutions of the resonator-qubit 
composite system. Similar method has been discussed 
in the context of q uantum computing with thermal ion- 
trap 48, 4^, 13, Hf and Josephson charge qubit [s^l- The 



feature of this coupling is the insensitivity to the quan- 
tum state of the resonator. In ion-trap quantum comput- 
ing, it is known as S0renson-M0lmer gate and has been 
experimentally demonstrated [H, Issl. ISSi? ]. However, 
the original S0renson-M0lmer relies on virtual excitations 
of the vibrational modes, whereas here the quantum bus 
carries real excitations. 

If a dc-pulsc is applied to to shift it from n(f>o, 
the time evolution of the composite system is driven by 
the Hamiltonian Eq. (|2ip . The operators included in the 
interaction Hamiltonian (a + at)^^'-* and the free Hamil- 
tonian (a+a^, iJk ■*) may be enlarged by their commutator 
into a closed Lie algebra of finite dimension. Thus the 
exact solution of the time evolution can be decomposed 
into a product over exponentials of the generators [56| . 
In the interaction picture, 



Hi 



(ate-* 



(23) 



The corresponding closed Lie algebra is 
^aax \ ax \ ai^^ (7^3^\ 1^ . The time evolution op- 
erator as the product of their exponentials can be 
written as 



C/i(t) 



^-'iD{t}^-iA(t)cr^icr^2 . 





24) 



This adds a 



9' 



a' 



term to the 
I with = 



Hamiltonian Eq. ([2T|l : 
-6M'^i^Ipy/uj/2L. How- 
ever since the qubit is far-detuned (e.g. according to the 
parameters used in Sec. V, A w 15.28 GHz and lo ^ 1 
GHz), this last term is a fast- rotating one and has neg- 
ligible contribution. In the following, we adopt Eq. ((2T|) 
as the effective system Hamiltonian. 



THE STRATEGY TO ACHIEVE EFFECTIVE 
TWO-QUBIT INTERACTION 

In this section, we discuss about how to achieve the 
two-qubit coupling in this composite system. The qubits 
only interact with each other indirectly through a com- 
mon quantum bus. In the dispersive limit, the operation 
time of two-qubit logic gate is limited by the small ra- 
tio g/5ijj where g is the qubit-bus coupling and 5uj the 
qubit-resonator detuning. In this case, the resonator is 
only virtually excited. In this paper we rely on another 
stradegy that the quantum bus carries real excitations. 
The effective two-qubit coupling is achieved by one or a 



where 




l)-t) 



(25) 



l)-t) 



In the following discussion, we neglect the universal phase 
factor D (t) . If the last factor of Eq. p4|) can be effec- 
tively canceled, 



Ui (t) = exp[-iD {t) - iA{t)axi<Jx2], 



(26) 



which represents the time evolution which is effectively 
governed by Hamiltonian ^ (TxiO'x2- This can be done in 
two different ways as described below: 



Single pulse operation 



By controlling the pulse length, a two-qubit gate is 
realized with a single dc-pulse which shifts <f>j |, from n$o 
(Fig. \5\ (a)). While the whole time evolution Eq. 
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is non-periodic, Eq. (|25|) shows that Bi{t) is a periodic 
function of time and it vanishes at r„ = 2mr/u!. At 
these times, the time evolution operator in the interaction 
picture reduces to 

C/i (T„) = cxp (i4n7r«;(i)K(2)^^^^^^^ (27) 

This is equivalent to a system of two coupled qubits with 
an interaction Hamiltonian cx (7xi<7x2- 

The minimum time to realize a rotation Uxx{0) = 
exp {i9axi(Tx2) is 



with 



mo = [no] 



2mQTT /ld 



6lj 



47rK(i)K(2) 



(28) 



(29) 



where [...] represents the integer part of a number. Note 
that we can not achieve a two-qubit rotation precisely un- 
less no happens to be an integer, so that no = mo- The 
error of one two-qubit gate is of the order of 47rK'^^)K'^^^ /lo 
(about 1% using practical parameters). This operation 
error can be avoided using a double-pulse method dis- 
cussed below. It is notable that increasing the frequency 
of the resonator lo cannot achieve a faster two-qubit gate 
(due to the y/ui dependence of g*-*-*), but it can reduce the 
error of the two qubit gate. 



Double pulse operation 

Alternatively a two-qubit logic operation can be con- 
structed with two successive operations as shown in 
Fig.EIb). 

Initially <I>j ^ is biased at n<I>o. The first dc-pulsc shifts 
it to a certain ^ n$o for a duration t/2. The evolu- 
tion operator (in the interaction picture) is Ui (|). After 
time t/2, one reverses the direction of the magnetic flux 
in the dc SQUID loop so that $j ^ is changed into -<I>d|, 
and f?''*-* = — g^'-*. The system is then driven by a new 
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FIG. 5: (Color online) Schematics of the pulse sequence to 
realize two-qubit gate operation. 



Hamiltonian H' = (— g^*)) for another t/2. The time 
evolution during the second pulse is 



U'j{t) = e-'^ We-'A (t)a.i<T.2 . 



n 



n 



with 



BS) = -B,it) 



A 



(30) 



(31) 



Note that the two terms in the brackets of Eq. (PD)) 
are permuted comparing with Eq. (|24p . so that the ex- 
pressions A'{t), D'{t) are different from A{t), D{t) in 
Eq. (ESI). 

The dynamics for the above two consecutive steps is 



C/tot it) 



(32) 



e-xp[-iM{t)axicrx2] exp[-z7V (t)] , (33) 



where 



M(t) 
Nit) 



25(1)3(2) /2 tut 
— sm t 

to \LU 2 

(g(l))2^(g(2))2 ^2 

— sm t 

LO 2 



(34) 



Therefore, the total time evolution is equivalent to 
the time evolution governed by two-qubit interaction 
~ cri^Vi^'' together with a universal phase factor. 

The time T to realize a rotation Uxxi^) in this way 
satisfies the nonlinear equation 



LoT 



uoT 



Olo 



2 4k(i)k(2)' 
In the case of a; ^ the solution is written as 

T ' 



2At(l)K(2) 



(35) 



(36) 



The two-qubit operation time is estimated using exper- 
imental parameters. Assuming Md — 10 pH and L = 100 
pH, one gets g''^^ ~ 36.02 MHz. As a cost of the low 
fluctuation related to the dc SQUID loop, the coupling 
strength associated with the dc SQUID loop is weaker 
than that with the main loop. For example, to realize a 
Uxx (f ), the operation time is about 204 ns. It is smaller 
than the qubit coherence time at the optimal point. The 
operation time is proportional to L/M^. Increasing the 
mutual inductance between the dc SQUID and the res- 
onator reduces the operation time. It is worth to point 
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out that the ratio g/ui is not required to be small. There- 
fore there is no fundamental limit on the operation time 
except the realizable coupling strength. 

As discussed in Sec. Ill, arbitrary single qubit gate 
can be performed after switching off the qubit-bus inter- 
action. Any non-trivial two-qubit gate can be built up 
with this XX coupling plus single qubit gates. For ex- 
ample, the C-phase gate can be constructed as (up to a 
global phase factor) [l| 



i?(e)^C/W(-^)L/p)(-^,cxp 



(37) 

with U'i\e) = exp (i6'o-i'V2) = c^^[iea'i^ /i^. Here 
we change representation so that = cTx and Gx ~ —Oz- 
And the CNOT gate can be readily constructed with C- 
phasc gate as 



CNOT = H^'^^Riir) H^'^\ 



(38) 



where H^"^") denote the Hadamard transformation on the 
i-th qubit as H^''^ = exp (^-in (ai'^ + ai'^) /2Vtj (up to 

a phase factor). 

With arbitrary single-qubit rotation and any non- 
trivial two-qubit rotation, universal quantum comput- 
ing can be realized according to quantum network theo- 
rem 1571. 



THE FEATURES OF THIS COUPLING 
PROTOCOL 



In the previous section, we have presented the way to 
realize two-qubit coupling and logic gate with our pro- 
posed setup. In this section, the features of this coupling 
protocol are analyzed with emphasis on the experimen- 
tal implementation. The qubit-qubit effective coupling 
commutes with the free Hamiltonian of the single qubit. 
This feature enables many practical advantages: 

(1) The main idea to implement a two-qubit operation 
from the exact evolution operator Eq. ([24|) is to cancel 
the part related with the degree of freedom of the res- 
onator, so that the final operation Eq. (|26p represents 
a qubit-qubit operation without entanglement with res- 
onator mode. Therefore the resonator mode does not 
transfer population with the qubit although the resonator 
mode mediates the qubit-qubit interaction. As a result, 
this two-qubit lo gic g ate is insensitive to the initial state 
of the resonator [49|. This feature is important for the 
experiment performed at finite temperature because the 
equilibrium state of the resonator is a mixed state. For 
example, there is 16% population at the excited state for 
a 1 GHz resonator at 30 mK. 

(2) As we mentioned, our coupling protocol works 
at the low-decohcrence optimal point where the qubit 
is robust to flux fluctuation and has long dccohercncc 



time. This is in contrast to other coupling protocols 
with dc-pulse control [H, li^. During the two-qubit 
operation, the control parameter is not the total mag- 
netic flux but rather a component in the dc SQUID loop. 
Therefore, the qubit can be biased at the optimal point 
<i>4*' = {n + l/2)4>o during two-qubit operation. 

While the dc SQUID adds a second control to the cir- 
cuit, it introduces extra decoherence. The fluctuation 
of the flux threading the dc SQUID loop results in the 
fluctuation of the energy splitting A and introduces de- 
coherence to the qubit dynamics. Suppose the magnetic 
flux are perturbed by the same amount of fluctuation 



as $ 



(i) 



(i) 



6^. Therefore the 



first-order effect of the fluctuation of magnetic flux in 
the main loop and the sub-loop are SE^ 



(i) _ JO 



6E'^'' = Cj'fJ^ respectively, with c^*'' = 



cl ' and 



and 



, where i?*-*' is the energy level spac- 



ig of the qubit, E^''> ^ Je"^ + A2 ($ 



If a 

qubit with E^/h = 259 GHz, Ej/E^ = 35 and 2ao = 
0.8 is biased at <i>t = $0/2 and 2ao cos (7r$d/*o) ~ 0.65, 
we get 



Si) 



163 GHz/$o 



(39) 



However, if qubit is not biased at the optimal point but 
close to the optimal point, e.g. e/E = 0.5, 



1100 GHz/$o 



(40) 



where we assume Ip = 500 nA. The influence of the fluc- 
tuation on the total magnetic flux is one order of mag- 
nitude larger than that on the dc SQUID loop. This 
suggests that although the dc SQUID loop introduces ad- 
ditional fluctuation to the system, the decoherence com- 
ing from flux fluctuation in dc SQUID is much less than 
that caused by shifting away from the degeneracy point 
<^f^ = $0/2. 

(3) A scalable qubit-qubit coupling scheme should al- 
low the coupling to be switched on-and-off (i.e. tun- 
able over several orders of magnitude). Otherwise, ad- 
ditional compensation pulse is needed to correct the er- 
ror in single-qubit operation. In our coupling protocol, 
as shown in Eq. (|19p . the external magnetic flux in the 
dc SQUID loop can be used to switch off the coupling by 
setting = 2n7r$o- When the qubit is decoupled from 
the data bus, single qubit operation can be controlled by 
^q^o independently. 

Our protocol does not require to change the amplitude 
of a dc pulse instantaneously. Finite rising and falling 
times of the controlling dc pulse will not induce addi- 
tional error to the two qubit coupling. This is essen- 
tially due to the qubit-resonator interaction commutes 
with the free Hamiltonian of the qubit at the optimal 
point. In the previous discussion, we assumed a constant 
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(7^''' for simplicity. In the experiments, the modulation 
of the magnetic flux always needs finite rising time, i.e.. 



g^^^' (t). As long as g^^' is a slow- varying (compar- 



ing with e~'^*) function of time t, the above discussion 
still holds except that the length of the pulse, i.e. T 
should satisfy 



B, (T) EE e-''^^.g(') (T) - g'-'^ (0) = 



(41) 



instead of T = 2mr/uj. The magnitude of the effective 
two-qubit interaction, i.e., A{t) in Eq. ((24|) is modified as 

AiT) = /'^^?i{e-*(g(i)Wg(2)(0)+5«(0)g(^)(i)) 

(42) 

To realize a certain xx rotation U = exp (^iOai^^ a^^^ is 

to apply a pulse satisfy Bi (T) = and A (T) ^ 9 simul- 
taneously. It is notable that the two conditions arc only 
related to the integral over the whole pulse and thus ro- 
bust to operation error. This conclusion is also applicable 
to the double-pulse method. 

(4) The evaluation is applicable to "ultra-strong cou- 
pling" regime where the coupling strength is even com- 
parable to the free Hamiltonian frequency as long as the 
approximation (jlSp is valid. Hence in principle, the two- 
qubit operation can be made as fast as single qubit op- 
eration. 

DISCUSSION AND CONCLUSIONS 

We illustrate two possible ways to scale up the two- 
qubit system. In Fig. [6] (a), a nearest-neighbor coupled 
qubits are sketched. They form a transverse Ising chain 
which can be used to implement quantum state trans- 
fer HI, 5^, 6^ and quantum information storage [6l| . It 





(a) 






(b) 



FIG. 6: (Color online) Schematics to scale up the coupling 
system, (a) Each qubit is coupled with the nearest neighbors 
by twisted LC resonators, (b) All qubits are interacting with 
a common TLR resonator on top of the qubits array. 



is possible to extend this configuration to 2D Ising model. 
Fig. [6] (b) shows an example to realize selective coupling 
between multiple qubits by a single quantum bus (such 
as a transmission line resonator). 

For nominally same parameters, there is natural spread 
of the junctions critical currents. This coupling mech- 
anism is robust to the difference of o'q^ and a^f^ be- 
cause the free Hamiltonian commutes with the interac- 
tion Hamiltonian. As such, in the sample fabrication pro- 
cess, the requirements on homogeneity and reproducibil- 
ity can be relaxed and meet with current production tech- 
nology. The additional on-site control lines require only 
one more layer. 

The qubit-resonator interaction commutes with the 
qubit free Hamiltonian. This feature enables quantum 
non-demolition (QND) measurement on superconducting 
qubit biased at the optimal point [g^I- This QND mea- 
surement is realizable even in the ultra-strong coupling 
limit [il. 
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